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Abstract. Many of the known extrasolar planets are "hot 
Jupiters," giant planets with orbital periods of just a few days. 
We use the observed distribution of hot Jupiters to constrain the 
location of the "inner edge" and planet migration theory. If we 
assume the location of the inner edge is proportional to the Roche 
limit, then we find that this edge is located near twice the Roche 
limit, as expected if the planets were circularized from a highly ec- 
centric orbit. If confirmed, this result would place significant limits 
on migration via slow inspiral. However, if we relax our assump- 
tion for the slope of the inner edge, then the current sample of hot 
Jupiters is not sufficient to provide a precise constraint on both the 
location and power law index of the inner edge. 



1. Introduction 

Early radial velocity discoveries were interpreted as showing a pile-up 
at an orbital period of three days, but recent transit surveys and very 
sensitive radial velocity observations have discovered planets with even 
shorter orbital periods. These discoveries suggest that the inner edge 
for hot Jupiters is not defined by an orbital period, but rather by a tidal 
limit which depends on both the semi-major axis and the planet-star 
mass ratio (See Fig. 1). This would arise naturally if the inner edge 
were related to the Roche limit, the critical distance at which a planet 
fills its Roche lobe. The Roche limit, a/j is is defined by 
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where Rp is the radius of the planet, Mp is the mass of the planet, and 
is the mass of the star. 

The numerous mechanisms proposed to explain the migration of 
giant planets to short period orbits can be divided into two broad cat- 
egories: 

1. Mechanisms involving slow inspiral, such as migration due to a 
gaseous disk or planetesimal scattering (Trilling et al. 1998, Gu 
et al. 2003). These would result in a limiting separation equal to 
the Roche limit. 

2. Mechanisms involving the circularization of highly eccentric orbits 
with small pericentcr distances, possibility due to planct-planct 
scattering (Rasio Sz Ford 1996, Ford, Havlickova, & Rasio 2001, 
Papaloizou Sz Terquem 2001, Marzari Sz Weidenschilling 2002), 
secular perturbations from a wide binary companion (Holman, 
Touma &; Tremaine 1997), or tidal-capture of free-floating planets 
(Gaudi 2003). These would result in a limiting separation of twice 
the Roche hmit (Faber et al. 2004). 

2. Statistical Analysis 

To quantitatively explore the observational constraints on the distribu- 
tion of hot Jupiters, we employ the techniques of Bayesian inference. In 
the Bayesian framework, the model parameters are treated as random 
variables which can be constrained by the actual observations. There- 
fore, to perform a Bayesian analysis it is necessary to specify both the 
likelihood (the probability of making a certain observation given a par- 
ticular set of model parameters) and the the prior (the a priori proba- 
bility distribution for the model parameters). Let us denote the model 
parameters by 6 and the actual observational data by d, so that the 
joint probability distribution for the observational data and the model 
parameters is given by 

p{d,e)=p{e)p{d\e)=p{d)p{e\d), (2) 

where we have expanded the joint probability distribution in two ways 
and both are expressed as the product of a marginalized probability dis- 
tribution and a conditional probability distribution. The prior is given 
by p{6) and the likelihood by p{d\9). On the far right hand side, p{d) is 
the a priori probability for observing the values actually measured and 
p{0\d) is the probability distribution of primary interest, the a poste- 
riori probability distribution for the model parameters conditioned on 
the actual observations, or simply the posterior. The probability of the 
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Figure 1.: Minimum mass ratio versus orbital period for the current 
observed sample of planetary companions. Planets discovered by radial 
velocity surveys are shown as triangles with arrows indicating 1 — a un- 
certainties in mass due to unknown inclination. The magenta squares 
have inclinations and radii measured via transits. The blue squares 
show planets discovered by transit searches. The green lines show 
the minimum mass corresponding to various velocity semi- amplitudes 
and roughly indicate where radial velocity surveys are nearly complete 
(> 30m/s), have significant sensitivity (> lOm/s), and are only begin- 
ning to detect planets (> 3m/s). The two red lines show the the location 
of the Roche limit (api) and the ideal circularization radius (fldrcj for a 
planet with a radius, Rp = 1.2Rj. The inner edge for the distribution 
of hot Jupiters is near Ocirc- Note that the red lines do not apply to 
the lowest mass planets that likely have a radius significantly less than 
\.2Rj due to their qualitatively different internal structure. 
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observations p{d) can be obtained by marginalizing over the joint prob- 
ability density and again expanding the joint density as the product of 
the prior and the likelihood. This leads to Bayes' theorem, the primary 
tool for Bayesian inference, 

^^^l'^^- p{d) - J dep{d\e)p{e) 

Often the model parameters contain one or more parameters of par- 
ticular interest (e.g., the location of the inner cutoff for hot Jupiters 
in our analysis) and other nuisance parameters which are necessary to 
adequately describe the observations (e.g., the fraction of stars with hot 
Jupiters in our analysis). Since Bayes' theorem provides a real proba- 
bility distribution for the model parameters, we can simply marginalize 
over the nuisance parameters to calculate a marginalized posterior prob- 
ability density, which will be the basis for making inferences about the 
location of the inner cutoff for hot Jupiters. 

Wc construct models for the distribution of hot Jupiters and use 
Bayes' theorem to calculate posterior probability distributions for model 
parameters given the orbital parameters measured for extrasolar planets 
discovered by radial velocity surveys. For the sake of clarity, we start 
by presenting a simplistic one-dimensional model for the distribution of 
hot Jupiters. We then gradually improve our model to understand how 
each model improvement affects our results. 

The primary question which we wish to address in this paper is the 
location of the inner edge of the distribution of hot Jupiters relative to 
the location of the Roche limit. Therefore, wc define x = a/aji, where 
a is the semi-major axis of the planet and qr is the Roche limit. We 
assume that the actual distribution of x for various hot Jupiters is given 
by a truncated power law, 

p{x\-f,xi,Xu)dx = x'^ (—\ , Xl < X < Xu, (4) 



X ) 

and zero else where. Here 7 is the power law index and xi and 
the lower and upper limits for x. The lower limit, x;, is the model 
parameter of primary interest, while 7 and Xu are nuisance parameters. 
Therefore, our results are summarized by the marginalized posterior 
probability distribution for x;. 

In order to minimize complexities related to the analysis of a pop- 
ulation, we choose to restrict our analysis to a subset of the known 
extrasolar planets for which radial velocity surveys are complete and ex- 
tremely unlikely to contain any false positives. To obtain such a sample, 
we impose two constraints: P < Pmax; where Pmax is the maximum or- 
bital period, and K > -fTmin, where -ftTmin is the minimum velocity semi- 
amplitude. We use Krain = 30m/s, based on the results of simulated 
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radial velocity surveys (Gumming 2004). We typically set Pmax =30d, 
even though radial velocity surveys are likely to be complete even for 
longer orbital periods (provided K > Kami)- This minimizes the chance 
of introducing biases due to survey incompleteness or possible structure 
in the observed distribution of planet orbital periods at larger periods. 
By considering only planets with orbital parameters such that radial ve- 
locity surveys are very nearly complete, our analysis does not depend on 
the velocities of stars for which no planets have been discovered. Note 
that our criteria for including a planet may introduce a bias depending 
on the actual mass-period distribution. We will address this point with 
a two-dimensional model at the end of our analysis. Also, our criteria 
exclude any planet discovered via techniques other than radial veloci- 
ties (e.g., transits), even if subsequent radial velocity observations were 
obtained to confirm the planet. 

Initially, we make several simplifying assumptions to make an an- 
alytic treatment possible. We assume uniform prior probability distri- 
butions for each of the model parameters, p{^) ~ f/(7min! Tmax) and 
p{xi,Xu) ~ const, provided xu < xi < Xu < Xuu and zero otherwise. 
The lower and upper limits (x;; and x^u ) for each parameter are chosen 
to be sufficiently far removed from regions of high likelihood that the 
limits do not affect the results. We assume that the orbital period (P), 
velocity semi-amplitude (K), semi-major axis (a), stellar mass (M*), 
and planet mass times the sin of the inclination of the orbit relative to 
the line of sight (m sin i) are known exactly based on the observations. 

We begin by assuming that sini = 1 for all planets and that all 
planets have the same radius, Rp. With these assumptions, the poste- 
rior probability distribution is given by 

n 

p{xi,Xu,j\xi,...Xn) ~ 7"(a;2 - a;/)"" n ^J~^ (^) 

provided that xu < xi < < < Xu < Xuu and 7min < 7 < 
7max- Here n is the number of planets included in the analysis, 
is the smallest value of x among the sample of hot Jupiters used in 
the analysis, and is the largest value of x in the sample. The 
normalization can be obtained by integrating over all allowed values of 
xi, Xu, and 7. 

We show the marginal posterior distributions in which we have 

integrated over the nuisance parameters, a;„ and 7 in Fig. 2 (left, dotted 
line), assuming Rp = 1.2Rj. The distribution has a sharp cutoff at Xf^ij 
and a tail to lower values reflecting the chance that xi < X(i-) due to the 
finite sample size. 

Next, we assume that the actual distribution of orbital inclinations 
is isotropic (cosi ~ C/[— 1,1]). For planets which were discovered by 
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Figure 2.: Left: Marginal posterior distribution for xi, the location of the 
inner edge for the population of hot Jupiters. The dotted line assumes 
all systems are viewed edge on (s'mi = 1). The solid line assumes an 
isotropic distribution of orbital inclinations ( except for transiting plan- 
ets). The remaining lines replace the assumption of R = l.2Rj for all 
non-transiting planets with a normal distribution for the planet radii 
using a dispersion dispersion of a^^p = 0.05i?,/ (long dashes), O.lRj 
(dots-dashes), or0.2Rj (short dashes). Upper Right: Dependence of the 
marginal posterior distribution for xi on the assumed mean planet ra- 
dius. {Rp) = l.ORj (long dashes), l.lRj (dotted), 1.2Rj (solid), l.SRj 
(dotted dashed), and 1.4i?j (short dashes), all assuming a^p = O.lRj. 
Lower Right: Same as above, but using a 2-d model (period & mass) 
which accounts for selection effects. 



radial velocities and the inclination was subsequently determined with 
the detection of transits, we use the measured inclination. The marginal 
posterior distribution for xi is shown in Fig. 2 (left, solid line). The 
sharp cutoff at a:;(i) is replaced with a more gradual tail, reflecting the 
chance that sini < 1 for planets with the smallest values of x. 

Next, we consider the consequences of allowing for a distribution 
of planetary radii. For transiting planets we use a normal distribution 
for the radius based on the published radius and uncertainty. For non- 
transitting planets, we assume a normal distribution of planetary radii 
with standard deviation, URp . We show the resulting marginalized pos- 
terior distributions in Fig. 2 (left). Allowing for a significant dispersion 
broadens the posterior distribution for xi and results in a slight shift to 
smaller values. 

We have also explored the effects of varying the model parameter 
-fmax) exploring values from 8d to 60d. We find that this does not make 
a discernible difference in the posterior distribution for xi. 
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Our results are sensitive to our choice for the mean radius for the 
non-transiting planets. In Fig. 2 (upper right) we show the posterior 
distributions for various mean radii, assuming cjRp = O.lRj. Since 
few planets have a known inclination, there is be a nearly perfect de- 
generacy between Rp and xi. Even when we include transiting plan- 
ets, this degeneracy remains near perfect, i.e., p{xi\Rp,xi, ...,Xn) — 

p{xi ■ -j^lR'pjXi, ...,Xn)- However, it can be seen that is extremely 
unlikely for xi to be near unity for any reasonable planetary radius. 

We have also performed an improved analysis using a two dimen- 
sional model which considers the joint planet mass-period distribution 
function. This allows us to account for observational selection biases 
due to the minimum mass for detecting a planet depending on the or- 
bital period. We assume that the distribution function is a truncated 
power law in both planet-star mass ratio and period. That is 



provided //mm < fJ, < /Umax, P < ^'max, and a{P, M^) > XI ■ aR{Rp,IJ,). 

Here fi = Mp/M^, and a and P are the new power law indices. We find 
the marginal posterior distribution for xi is very similar to the results 
of our 1-d analysis. The most significant diff'erence is that the posterior 
distribution for xi shifts slightly to wards larger separations (see Fig. 2, 
lower right). 

3. Discussion 

The current distribution of hot Jupiters discovered by radial velocity 
searches shows a cutoff that is a function of orbital period and planet 
mass. Our Bayesian analysis solidly rejects the hypothesis that the 
cutoff occurs inside or at the Roche limit, in contrast to what would be 
expected if the hot Jupiters had slowly migrated inwards on a circular 
orbit. Instead, our analysis shows that this cutoff occurs at a distance 
nearly twice that of the Roche limit, as expected if the hot Jupiters 
were circularized from a highly eccentric orbit. These findings suggest 
that hot Jupiters may have formed via planet-planet scattering (e.g., 
Rasio &; Ford 1996), tidal capture of free floating planets (Gaudi 2003), 
or secular perturbations from a highly inclined binary companion (e.g., 
Holman, Touma & Trcmaine 1997). If the hot Jupiters indeed were 
circularized from a high eccentricity orbit, then this raises the challenge 
of explaining the origin of giant planets with orbital periods ~ 10 — 100 
days. 

An alternative explanation is that the planets migrated inwards on 
a nearly circular orbit at a time when the planets were roughly twice 
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their current radii. Future observations of low mass planets may make 
it possible to test this alternative, assuming that the time evolution of 
their contraction is significantly different than for Jupiter-mass planets. 
A third alternative is that short-period giant planets are destroyed by 
another process before they reach the Roche limit. HST observations of 
HD 209458 indicate absorption by matter presently beyond the Roche 
lobe of the planet and have been interpreted as evidence for a wind 
leaving the planet powered by stellar radiation (Vidal-Madjar ct al. 
2003, 2004). Further theoretical work will help determine under what 
conditions these processes can cause significant mass loss. 

Future planet discoveries will either tighten the constraints on the 
model parameters or provide evidence for the existence of planets defi- 
nitely closer than twice the Roche limit. In particular, new discoveries 
of very low mass planets could better constrain the shape of the inner 
cutoff as a function of mass. In the future, an improved analysis could 
also include such low-mass planets where surveys are not yet complete. 
For future theoretical work, we hope to explore the possibility of orbital 
circularization occurring at larger orbital radii, possibly in a protoplan- 
etary disk or while the star is young and still contracting. 
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